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Abstract 

With nonrelativistic QCD factorization, we calculate the T(IS') prompt production at hadron 
colliders at next-to-leading order in a s . In addition to the color-singlet contribution, color-octet 
(CO) channels (especially the P-wave channel) up to 0(v ) are all considered. Aside from direct 



a; 

production, the feed-down contributions from higher excited S-wave and P-wave bb states to T(15) 



production are also included. We use the potential model estimates as input for color-singlet 
(CS) long-distance matrix elements (LDMEs). While for CO contributions, we find they can be 
approximately described by three LDMEs: (0( 3 sf ] )), (0(!s{) 8] )) and (ofpj 81 )). By fitting the 
Tevatron data we can determine some linear combinations of these LDMEs, and then use them to 
predict T(1S) production at the LHC. Our predictions are consistent with the new experimental 
data of CMS and LHCb. 

PACS numbers: 12.38.Bx, 13.85.Ni, 14.40.Pq 
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I. INTRODUCTION 



The study of heavy quarkonium production is particularly interesting because it may 
provide decisive information in understanding hadronization of heavy quarks and gluons in 
QCD. The most widely accepted theory to describe heavy quarkonium production at present 
is nonrelativistic QCD (NRQCD) factorization [I], in which the production is factorized 
into perturbative calculable short- distance coefficients and nonperturbative (and universal) 
long-distance matrix elements (LDMEs). As short- distance coefficients can be expanded in 
strong-coupling a s and each LDME has a definite power in v (the velocity of heavy quarks 
in the rest frame of heavy quarknonium) , NRQCD factorization gives predictions by dou- 
ble expansion in a s and v 2 . NRQCD factorization is efficient such that in principle only 
a finite number of universal parameters, which can be determined by using some known 
experimental data, are involved with required precision in predicting other production pro- 
cesses . Although a complete proof of factorization is still lacking, at least it holds up to 
next-to-next-to-leading order (NNLO) in a s 0,0. 

Based on NRQCD factorization, charmonium production in hadron colliders has been 
studied extensively in recent years PHI?]. Specifically, for J /if) hadroproduction, it is found 
that all data at large pt, including both yield and polarization, can be well described by 



NRQCD factorization if one chooses a large M and a small Mi [Tj| (see also refs.[[12|, |Hfl). 
Here M and Mi are linear combinations of related LDMEs which are defined in refs. fT2|, 
I4| and will also be mentioned below, and roughly speaking, their values represent the 



importance of p^ 6 behavior and p^, 4 behavior in J / ip production cross sections, respectively 
There are reasons that studying bottomonium may be a more suitable choice than char- 
monium to test the NRQCD factorization formalism. First, the value of v 2 is smaller in 
bottomonium (« 0.1) than that in charmonium (« 0.3), thus the expansion in v 2 should 
converge faster in bottomonium. Second, the mass of bottomonium is about 3 times of that 
of charmonium, then asymptotic freedom implies the convergence of a s expansion is also 
better in bottomonium. However, on the experiment side, the situation is not so satisfactory 
as the production rates of bottomonium are much smaller than charmonium, e.g., the cross 
section of T(1S) is about two orders of magnitude smaller than that of J/if>. Furthermore, 
there are more excited bottomonium states which are below the open bottom (say BB) 
threshold and can decay into lower bottomonium states such as T(IS') with large branching 
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ratios and consequently contribute a substantial fraction to the lower bottomonium inclusive 
production by the so-called feeddown contributions, thus it is hard to measure the direct 
production from the prompt inclusive production. In the LHC era, we expect these disad- 
vantages may be overcome by the higher luminosity, thus testing NRQCD factorization by 
bottomonium production seems to be hopeful. 

The inclusive differential cross section and polarization of T are measured at the Teva- 
tron [p2|-p5|], but the T(1S) polarizations observed by DO f24| and CDF [|25| disagree with 
each other. Furthermore, both the DO and CDF measurements contradict the LO NRQCD 
prediction p6 |. As argued in refs. [iPT[ [14], the NNLO and even higher-order contributions 
pTfl may not be important, as compared with the full NLO QCD contributions including 
both color-singlet (CS) and color-octet (CO) channels, which are essential in understand- 
ing the T (and similarly the J/ip) hadroproduction. Partial NLO QCD contributions to T 
hadroproduction have been calculated recently ||, [7|, ||, |28|, and it is found that the NLO 
QCD corrections of S-wave CO channels only slightly change the transverse momentum dis- 
tribution and the polarization, while the correction of CS channel may bring on significant 
enhancement to the momentum distribution and change the polarization from transverse 
at LO into longitudinal at NLO. But the NLO contributions of P-wave channels for T 
hadroproduction are still missing. 

At the LHC, CMS has published the first run data for T production P9| , and LHCb has 
also reported the measured result |JIJ. Thus it is timely to present a complete NLO theo- 
retical prediction for T production, and compare theory with experiment. In this work, we 
study the T(1S) hadroproduction in the framework of NRQCD, including all NLO contribu- 
tions and feed-down contributions. The paper is organized as follows. We briefly introduce 
our calculation in Sec. |J. Then in Sec. [II, we describe our method for taking into account 
the feed-down contributions. (Note that the feed-down contributions for T(1S) have not 
been treated seriously in all previous theoretical works.) In Sec. [IV], we fit data to determine 
LDMEs and then give predictions for the LHC experiment. Finally, a summary is given in 
Sec. 0. 
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II. NLO CALCULATION 



The method of NLO calculation used in the this work is similar to that used in J /if) and 



Xc production |TT], [12], |14|| . For completeness, we will sketch it in this section. 

According to the NRQCD factorization formalism, the inclusive cross section for direct 
bottomonium H production in hadron-hadron collisions is expressed as 



da\pp ^H + X] = J2 da[{bb) 



J2 / dx 1 dx 2 G i/p G j/p x da[i + j ->• (bb) n + X]{0%), 



where p is either a proton or an antiproton, the indices i, j run over all the partonic species, 
and n denote the color, spin and angular momentum (L n ) of the intermediate bb states. In 
this work, we calculate the cross sections up to v 4 corrections, so that the intermediate states 
include 3 Sj 1] , 3 P l j\ !sj 8] , 3 Sf ] and 3 pf\ Note that our definition of CS LDMEs (0 H { 3 Sf ] )) and 
(O h ( 3 P^)) are different from that in Ref. |l| by a factor of l/(2iV c ). The calculation proceeds 
with three steps: calculating the parton level differential cross section da[i + j — > (bb) n +X], 
integrating over the phase space, and fitting the LDMEs. 

NLO corrections for the parton level differential cross section include virtual corrections 
and real corrections. For virtual corrections, we use FeynArts (3l| to generate Feynman 
diagrams and amplitudes. We then calculate these thousands Feynman diagrams analyti- 
cally using our self-written Mathematica code. Finally, we output the simplified expression 
into C++ code. Because the infrared (IR) divergence will appear when doing phase space 



integration for the real correction, we use the two cutoff phase space slicing method |]32] to 
isolate the divergence. The contributions from the singular phase space part are calculated 
analytically, while finite parts are calculated by using the Berends-Giele off-shell recursive 
relations P5| . 

In the analytical calculation we have checked that all the ultraviolet (UV) and infrared 
(IR) singularities are canceled exactly. The UV divergences are removed by renormalization. 
The IR singularities arising from loop integration and phase space integration of the real 
correction partially cancel each other. The remaining IR singularities are absorbed into the 
proton parton-distribution-functions (PDFs) and the NRQCD LDMEs. 

The numerical integration over the phase space is handled by our self-written C++ codes, 



where we also use both QCDLoop p4| and LoopTools [3q] to calculate the scalar functions 



in the virtual corrections numerically. We verified that our results are independent of the 
two cuts introduced by the phase space slicing method. The method of fitting LDMEs will 
be discussed in Sec. [TV[ 

III. TREATMENT OF FEED-DOWN CONTRIBUTION 

One difficulty in predicting T production cross section is the treatment of feed-down 
contribution. There are several higher excited states that can decay into T(1S) and they 
include: T(2S), T(3S), Xbi(lP), x 62 (l-P)> Xbi(2P) and Xb2^P)- Unfortunately, there are 
not enough data to determine LDMEs of these higher excited states, therefore, it is hard 
to predict their feed-down contributions to T(IS'). In fact, all previous predictions that 
were based on NRQCD factorization did not have a serious treatment of the feed-down 
contributions. 

The key point to deal with feed-down contribution is to determine the relation between 
momentum of higher excited states and momentum of T(IS'). In ref. [0, we find a very 
good approximation that the ratio of two momenta is inversely proportional to the ratio of 
their masses. Notice that the mass differences between these excited states and T(IS') are 
of the order of m^v 2 and v 2 is very small in bottomonium, as a result, unlike the J /if) case, 
the momentum shift can be ignored when these excited states decay into T(IS'). Hence the 
production LDMEs can be approximately combined into 6 independent ones: 

(£>( 3 Si 1] )> =(C T(15) ( 3 5! 1] )) + (O nnS) Cs [ i ] )}Bi(T(nS) ->■ T(1S)), 

n=2,3 

(0( 3 Pl 1] )) = (C X6l( " P) ( 3 P 1 [11 ))Br(x 61 (nP) T(15)), 

n=l,2 

(OCPP)) = (0^ nP \ 3 P?))BT( Xb2 (nP) -> T(1S)), 

n=l,2 

(0Csf ] )} =(C T(15) ( 3 M 81 )) + Yl (O r{nS) Cs [ i ] )}Bi(T(nS) -> T(1S)) (2 ) 

n=2,3 

+ E E^ Xftj(nP) W ] )) Br (^( np )^ T ( 15 ))' 

n=l,2 J=l,2 

(C(!si 8] )) =(C T(15) (!si 81 )) + (C ,T(n5) ( 1 5'(f 1 ))Br(T(nS) -> T(1S)), 

n=2,3 

(OCP^)} =(C T ( 15 )( 3 P ( i 81 )) + (C T(n5) ( 3 P(! 8] ))Br(T(nS) -)• T(1S)). 

n=2,3 
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Here the Xba(lP, 2-P) feeddown into T(IS') is ignored due to the smallness of the transition 
branching ratios. The potential model results of wave functions and their derivatives at the 
origin can be chosen as[|36| 



-Rr(15)(C 


)| 2 


= 6.477 GeV 3 


-Rt(25)(C 


)| 2 


= 3.234 GeV 3 


-Rr(3S)(C 


)| 2 


= 2.474 GeV 3 




)| 2 


= 1.417 GeV 5 


# Xf) (2P)( C 


)| 2 


= 1.653 GeV 5 



Air 
3 



\R 



T(nS) 



(3) 



(4) 



and the CS LDMEs can be estimated by 

(0^ (nP) #f)> = ^l^(nP)(0)| 2 (2J+ 1). 

With the PDG data of branching ratios [37], we get 

(C( 3 4 1] )) =1.81 GeV 3 
(0( 3 PP)) =0.54 GeV 5 
(OfpJ 11 )) =0.62 GeV 5 . 

Now there leave only 3 unknown CO LDMEs: (0{ 3 Sf ] )), (oCs 1 ^)), and (0( 3 P®)). They 
will be determined by fitting the Tevatron data fl22| , |23| . Because they are nearly universal 
(up to a correction of order v 2 with calculated short- distance coefficients in the fit), the 
fitted results can be used to predict T(1S) production in other colliders. 
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IV. NUMERICAL RESULT 

The CTEQ6L1 and CTEQ6M PDFs (parton-distribution functions) [j38|1 are used for LO 
and NLO calculations respectively. The bottom quark mass is set to be mj, = 4.75 GeV, 
while the renormalization, factorization, and NRQCD scales are fi r = fif = tjit and /ia = 
where my = \/ p\ + * s the T transverse mass. The center-of-mass energies are 1.8 TeV, 
1.96 TeV and 7 TeV for the Tevatron RUN I, RUN II, and LHC respectively. 

The Tevatron data [fZ2]. in Fig. [1] with pj- > 8 GeV are used to fit the 3 unknown 
CO LDMEs: (0( 3 Sf ] )), (O^Sf )) and (0^)). The reason for discarding the low p T data 
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FIG. 1: Transverse momentum distributions of prompt T(IS') production cross sections at the 



Tevatron. The CDF data are taken from Ref. [22]. The DO data are taken from Ref. 



is that these data are far from the large pt region <C 1), and may also be affected 



by nonperturbative effects, which can not be described by our fixed order perturbative 



calculation. As discussed in ref. |12|, |14| , by fitting large px data at the Tevatron one can 
only constrain two linear combinations of LDMEs that have p^ 4 and pT^ behaviors at parton 
level respectively. Thus, to have a constrained fit, the short- distance coefficient of P-wave 
channel is decomposed into linear combination of that of two S-wave channels in ref. ||1 2|| , 
and as a result, one needs to fit only two linear combinations of LDMEs (Mo and Mi). 
However, for T production we find the decomposition of P-wave channel is good only for 
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Pt > 15 GeV, thus, to fit the Tevatron data with px > 8 GeV, we cannot decompose the 
P-wave channel. Therefore, we fit the three LDMEs using a similar method described in 
ref . \M . 



Define 



and the correlation matrix C 



By minimizing \ 2 we have 



C 



2 = (oCsf ] )), 

3 = 



(O(W) 



mt 



( 



\ 



r -i_i rfV 

ij 2dO i dO j 



0.24 -0.024 -0.54 
-0.024 0.0025 0.0054 
-0.54 0.054 1.21 



The eigenvalues Aj with corresponding eigenvectors Vj of C are 



(6) 



(7) 



(8) 



Ai = 1.5, "v\ = (-0.41, 0.040, 0.91) 
A 2 = 3.5 x 10" 4 , "v^ = (0.79, -0.48, 0.38) 
A 3 = 1.3 x 10" 5 , ^ 3 = (0.46,0.87,0.17). 



(9) 



The LDMEs corresponding to the eigenvectors are 



/ 0l \ 



where we denote matrix 



Ai 
A 2 
V A 3/ 



V 



o 2 
v o 3/ 



V 



(10) 



(11) 
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Inserting Eqs.(§) and @ into Eq.(|T0"l), we have 

Ai = -274 x 1(T 2 GeV 3 (±44%), 

A 2 = 6.04 x 10~ 2 GeV 3 (±31%), (12) 
A 3 = 10.5 x 10~ 2 GeV 3 (±3.4%). 

In this way, the three CO LDMEs are expressed in terms of their linear combinations Aj, 
which correspond to the eigenvectors of the correlation matrix. Because the Tevatron data 
are not enough to determine all the three CO LDMEs, Ai is unphysically determined, which 
explains why the absolute value of Ai is much larger than A 2 and A3. Nevertheless, to 
change A x to be the same order as A 2 and A 3 , one needs more than two a shift. It implies 
results in Eq.(^) may underestimate the error of Ax. 

Values of Aj contain main result in our fit. To use them to predict T(1S) production in 
other experiment, we express the differential cross section as 

3 3 

da = ^^dfriOi = ^^ajAj, with "2" = da V' 1 , (13) 

i=l i=l 

where d&i denote corresponding short-distance coefficients. In this form, the errors induced 
by Aj can be easily taken into consideration for they are independent. Based on Eq. ([T3"D , 
our predictions for CMS and LHCb are plotted in Fig. ^| and Fig. |3|, respectively, where 
CMS and LHCb data are taken from refs. [^, |31J. The uncertainties of the curves con- 
cern the renormalization scale dependence in the calculation and the errors from A,. We 
treat these two types of uncertainties as independent ones. From these figures, we can see 
that our predictions are consistent with the LHC experimental data, which is an explicit 
demonstration of the universality of LDMEs defined in Eq. (||). 

V. SUMMARY 

In summary, we calculate the complete NLO corrections for the Y(lS) production at 
hadron colliders up to 0{a i s v i ). Ignoring corrections of higher-orders in v 2 , we combine the 
production LDMEs of T(IS') and other excited states into 3 color-singlet LDMEs and 3 
color-octet LDMEs. These 6 LDMEs are approximately universal and they include almost 
all feed-down contributions to T(IS') production. The CS LDMEs are estimated by using 
potential model results, while the CO LDMEs are determined by fitting the Tevatron data. 
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FIG. 2: Transverse momentum distributions of prompt T(1S) production cross sections at the 
LHC. The CMS data are taken from Ref.gjJ. 1Q 
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FIG. 3: Transverse momentum distributions of prompt T(1S) production cross sections at the 
LHC. The LHCb data are taken from Ref . J3Q[| .-\ -\ 



Then we find our predictions well coincide with the new experimental data at the LHC Our 
work may provide a new test for the universality of LDMEs in T(1S) hadroproduction. 

To have a comprehensive understanding of T(IS') hadroproduction, it is certainly impor- 
tant to also compare the theoretical result with the polarization data for T(1S), we leave it 
as a further study Encouraged by the result of J / ip polarization |T^] , where we find the J / ip 
polarization and yield can be consistently explained by two well constrained CO LDMEs 
(M and Mi), a good description for the T(1S) data including yield and polarization seems 
to be promising. However, note that the values of two well constrained CO LDMEs in 



J/tp production are significantly different [[12] , [L4J], while for T(IS') production we find from 



Eq.(|l2]) that A 2 and A 3 are of the same order. Another complexity concerns the influence of 
big feedown contributions on T(IS') polarization. Therefore, a full understanding of Y(IS') 
production including both yield and polarization may provide important information in ad- 
dition to the study of J /ip production. On the experiment side, because the bottom quark 
is heavy: mj m 5 GeV, to test the large pr (— "C 1) behavior one needs to measure the 
cross sections and polarizations at pt as large as, say 30 GeV and even larger, with higher 
statistics, and to separate the higher excited bb production from the Y(IS') production. This 
is a hard task for experiment, and we hope it can be fulfilled at the LHC in the near future. 
Then we can make more thorough comparison between theory and experiment, and provide 
a further test of NRQCD factorization. 
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